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Building materials are often incorporated 
into complex, multilayer macro structures 
that are simply not amenable to measure- 
ments using coax or waveguide sample 
holders. In response to this, we developed 
an ultra-wideband (UWB) free-field meas- 
urement system. This measurement system 
uses a ground-plane-based system and two 
TEM half-horn antennas to transmit and 
receive the RF signal. The material sam- 
ples are placed between the antennas, and 
reflection and transmission measurements 
made. Digital signal processing techniques 
are then applied to minimize environmen- 
tal and systematic effects. The processed 
data are compared to a plane-wave model 



to extract the material properties with opti- 
mization software based on genetic algo- 
rithms. 
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1. Introduction 

The NIST Time-Domain Fields Project has devel- 
oped a nondestructive, ultra-wideband (UWB) meas- 
urement system to measure complex permittivities of 
building materials. The purpose is to develop a data- 
base of material permittivities for use by engineers for 
modeling wireless signal propagation in buildings and 
related applications. 

The NIST system consists of an aluminum ground 
plane, two TEM horn antennas, precision microwave 
cables, and a vector network analyzer. A wall of sample 
building material is placed between the antennas, and 
both transmission and reflectivity are measured over a 
wide frequency range. The stepped-frequency S- 
parameters are then post-processed to isolate the sam- 
ple response and reduce systematic effects due to the 
environment. The net result is a set of gated S-parame- 
ters. Material properties are then extracted by an opti- 
mization process based on plane-wave assumptions and 



the gated S-parameters. This optimization process uti- 
lizes a genetic algorithm (GA) [1] that provides accu- 
rate estimates of the electrical properties of a sample 
under test. This system has the marked advantage over 
previous systems in that it permits the measurement of 
very heavy wall samples without the need of a sample 
holder. The system also has high bandwidth, range 
(spatial) resolution, and measurement fidelity. 

This paper consists of six sections: Section 2 
describes the measurement system. Section 3 intro- 
duces the measurement sequence and the initial data 
processing. Section 4 describes the mathematical 
model and how the GA is used to optimize material 
parameter estimates for two candidate materials, cross- 
linked polystyrene and concrete. Section 5 compares 
the results obtained thus far with other measurement 
techniques and published results. Section 6 presents 
conclusions. 
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2. The Measurement System 

Figure 1 is a diagram of the laboratory-based meas- 
urement system. It consists of an aluminum ground- 
plane, two TEM half-horn antennas (ground-plane 
mounted), several material samples, two precision 
microwave interconnecting cables, and a two-port 
Vector Network Analyzer (VNA). Four S-parameters 
(S n , S 2 i, S 12 , S 22 ) are measured at 1601 stepped frequen- 
cy points in the range of 30 MHz to 6 GHz. For meas- 
urements of plywood, gypsum, and cross-linked poly- 
styrene, the 7.3 m x 7.3 m aluminum ground-plane sys- 
tem of the NIST cone and ground-plane time-domain 
range was used. Measurements of the concrete slab 
were made on a smaller 3.66 m x 3.66 m aluminum 
ground-plane. Sample thicknesses were measured with 
precision calipers at various points around the sample 
and are accurate to within ±0.1 cm. Figure 2 is a photo- 
graph of the measurement system being used on a sam- 
ple slab of cross-linked polystyrene, and Fig. 3 is a pic- 
ture of the concrete slab measurement setup. The TEM 
half-horns can be seen with their styrofoam wedge 
mounts on either side of the concrete slab. 



3. Measurement Sequence and Data 
Processing 

The measured S-parameters must be post-processed 
further before material properties can be extracted. In 



order to quantify the scattering from the sample under 
test, three measurements must be made: (1) a back- 
ground measurement with the two boresighted anten- 
nas at a specified separation, (2) a reference reflectivi- 
ty measurement with a large aluminum sheet halfway 
between the antennas, and finally, (3) the building 
material sample placed halfway between the antennas. 
The processing consists of three principal steps: (1) 
The background S-parameters, measured as a function 
of frequency with no sample in place, are subtracted 
from the S-parameters measured with either a sample 
or metal sheet in place; (2) these subtracted S-parame- 
ter sets are then inverse Fourier transformed, time 
gated, and Fourier transformed by NIST-developed 
software to yield a set of "gated" S-parameters; and (3) 
the gated S-parameters are then normalized with either 
the metal reference for reflectivity or with the antenna- 
to-antenna reference for transmission [2]. Specifically 
we have 
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Fig. 1. Plan view and profile view of sample under measurement. 
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Fig. 2. 1.2 m x 1.2 m cross-linked polystyrene in styrofoam holders for measurement on a laboratory 
ground-plane. 




Fig. 3. 1 .2 m x 1 .2 m x 0.2 m concrete measurement between TEM half-horns on portable ground-plane. 
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These three steps yield a set of normalized free-field S- 
parameters (5 n , S 2 i, S 12 , S 22 ) that quantify both the 
reflection and transmission through the sample under 
test. 



4. Extraction of Material Electrical 
Properties 

The free-field S-parameters are used in the determi- 
nant of the S-matrix [3]: 

F 2 _ 2 
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where 
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domly generating many possible solutions. These solu- 
tions are tested with the objective function and ranked 
by a minimum cost criterion. A certain fraction of the 
best solutions is selected and a new set of solutions is 
generated. This process is repeated until a desired level 
of convergence is met or a maximum number of itera- 
tions is achieved. A unique property of the GA is that 
per iteration, a fraction of the solutions are randomly 
changed to ensure that the results are not subject to 
local minima. The cost function that is used is the 
absolute distance per frequency between calculated and 
measured results and is given by 



S(£ r ,£ r )\ Mcalc) -\S 
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(9) 



The first term in Eq. (9) is equivalent to the right-hand 
side of Eq. (4), and the second term to the left-hand side 
of Eq. (4). The cost function of Eq. (9) is minimized. 



where /is the operating frequency, 
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In Eqs. (4) through (8), L is the measured thickness of 
the sample, d is the antenna separation without a sam- 
ple present, and c vac and c lab [4] represent the speed of 
light in vacuum and in the laboratory, respectively. 
Equation (8) contains the complex relative permittivity 
e* = s' r -js'r and complex relative permeability 
jU* = jx r -jii r of the sample. We assume all sample 
materials are nonmagnetic and do not have conductive 
properties, so that jx' r = 1 and [i° r = 0. 

Complex relative permittivity is determined by 
adjusting the value of e* at each frequency to optimize 
the least-squares difference between the measured 
results and the theoretical results predicted by Eqs. (4) 
through (8). The optimizer used for this purpose is a 
genetic algorithm (GA). In short, the GA [1] is a pow- 
erful technique used to minimize the difference 
between measured and calculated results and, from this, 
to extract the desired parameters. A GA works by ran- 



5. Results 

Complex permittivity results are presented for cross- 
linked polystyrene and concrete. In the case of cross- 
linked polystyrene, we measured a 1.2 m x 1.2 m sam- 
ple that was 2.9 cm thick. In this case we know, based 
on published data obtained by a Split-Post Resonator 
technique [5], that £ r ' = 2.55 at two frequency points, 
1.44 GHz and 2.06 GHz. Using this as a guide, the ini- 
tial population for the genetic algorithm of e r is drawn 
from the domain 1 < e r <3 and the results are shown in 
Fig. 4. This technique does not have sufficient sensitiv- 
ity to measure £ r for a low-loss material like cross- 
linked polystyrene. 

A series of measurements was also made on a 
1.2 m x 1.2 m x 0.2 m section of concrete wall. This 
sample weighs approximately 907. 1 8 kg and is an ideal 
candidate for ground-plane based measurements. The 
results of the concrete measurements are shown in Figs. 
5 and 6. The results are in good agreement with inde- 
pendent measurements, using a shielded open-circuit 
technique conducted by the Electrical Properties of 
Materials Project, and with published data [2]. The 
major sources of uncertainty are time-gating errors and 
possible edge effects that are not gated out of the meas- 
urement. Individual source uncertainties are presented 
in Table 1 where 8R and 8T are the standard deviations 
in reflection and transmission signal levels for the stat- 
ed uncertainties at a frequency of 1 GHz. These indi- 
vidual uncertainties, at each frequency, are RSS com- 
bined into reflection and transmission uncertainties, AR 
and AT. These uncertainties are propagated into equa- 
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Fig. 4. ef. versus frequency for a Split-Post Resonator method [4] and 1 .2 m x 1 .2 
m x 2.9 cm cross-linked polystyrene. 
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Fig. 5. e/ versus frequency for a Shielded Open-circuited method and 1.2 m x 1.2 
m x 0.2 m concrete. 
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Fig. 6. £ r " versus frequency for a Shielded Open-circuited method and 1.2 m x 1.2 
m x 0.2 m concrete. 



Table 1. Sources of uncertainty and their associated values at 1 GHz 
(*=1) 



Source of uncertainty 



Combined uncertainty 
(AR and AT) 



SR 



0.0381 



8T 



Repeatability 


0.0176 


0.0005 


Time Gating 


0.0187 


0.0213 


Edge Effects 


0.0274 


0.0045 


Length 


0.0049 


0.0049 


Phase Measurement 


0.0042 


0.0042 



0.0225 



6. Conclusions 

Utilizing a NIST-developed free-field measurement 
system, we were able to effectively extract electromag- 
netic parameters of materials with appropriate signal 
processing and post-processing using a genetic algo- 
rithm. In this process, an objective function chosen to 
represent the physics of the system is minimized, giv- 
ing very accurate results. The data obtained should be 
of great value to the UWB community, and with further 
research this system promises to be even more accu- 
rate, provided that the sources of uncertainty are 
addressed. 



tions for the permittivity uncertainties for reflection and 
transmission. These permittivity uncertainties are com- 
bined in an RSS uncertainty with a coverage factor of 1 
(k = 1). These uncertainties are shown by the error bars 
on each plot. All in all, the concrete results are very 
encouraging. Note also that concrete is highly disper- 
sive and lossy as understood by the imaginary part of 
permittivity — a result that should be of great interest to 
the wireless community. 



Acknowledgments 

The authors thank OLES for initially funding this 
research and Dr. Dennis Friday and Dr. Perry Wilson of 
the NIST Electromagnetics Division for continued sup- 
port. We also thank Dr. Edward Kuester of the 
University of Colorado at Boulder for much guidance. 



72 



Volume 112, Number 1, January-February 2007 

Journal of Research of the National Institute of Standards and Technology 



7. References 

[1] R. L. Haupt and S. E. Haupt, Practical Genetic Algorithms, Wiley 
Interscience, 2nd Ed., John Wiley & Sons, Inc. (2004). 

[2] R. T. Johnk, D. Novotony, C. A. Grosvenor, N. Canales, J. A. 
Graham, D. Martin, S. Yencer, D. Hibbard, L. L. Nagy, and T. 
Roach, An Electric Field Uniformity Study of an Outdoor 
Vehicular Antenna Test Range. Natl. Inst. Stand. Technol. Tech. 
Note, to be published. 

[3] J. Baker-Jarvis, M. D. Janezic, B. F. Riddle, R. T. Johnk, P. 
Kabos, C. L. Holloway, R. G Geyer, and C. A. Grosvenor, 
Measuring the Permittivity and Permeability of Lossy Materials: 
Solids Liquids, Building Material and Negative-Index Materials. 
Natl. Inst. Stand. Technol. Tech. Note 1536 (February 2005). 

[4] H. Liebe, An Updated Model for Millimeter Wave Propagation in 
Moist Air. Radio Science, vol. 20, pp. 1069-1089 (Sep.-Oct. 
1985). 

[5] NIST Standard Reference Material 2870, Certificate of Analysis. 

[6] R. T. Johnk, A. R. Ondrejka, and C. L. Holloway, Time-Domain 
Free-Space Evaluations of Urethane Slabs with FT-TD 
Computer Simulations. IEEE Int. Symp. Digest on 
Electromagnetic Compatibility (1998) pp 290-295. 

About the authors: Ben Davis, Chriss Grosvenor, 
Robert Johnk, David Novotny, James Baker-Jarvis, and 
Michael Janezic are electronics engineers in the 
Electromagnetics Division of the NIST Boulder 
Laboratories. The National Institute of Standards and 
Technology is an agency of the Technology 
Administration, U.S. Department of Commerce. 



73 



